The effect of electrostatic shielding using invisibility cloak 
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The effect of electrostatic shielding for a spherical invisibility cloak with arbitrary charges inside 
is investigated. Our result reveals that the charge inside the cloak is a crucial factor to determine 
the detection. When charged bodies are placed inside the cloak with an arbitrary distribution, the 
electric fields outside are purely determined by the total charges just as the fields of a point charge at 
the center of the cloak. As the total charges reduce to zero, the bodies can not be detected. On the 
other hand, if the total charges are nonzero, the electrostatic potential inside an ideal cloak tends to 
infinity. For unideal cloaks, this embarrassment is overcome, while they still have good behaviors of 
shielding. In addition, the potential across the inner surface of an ideal cloak is discontinuous due 
to the infinite polarization of the dielectric, however it can be alternatively interpreted as the dual 
Meissner effect of a dual superconductive layer with a surface magnetic current. 

PACS numbers: 41.20.Cv, 42.79.-e 



I. INTRODUCTION 

Since the pioneering works of Pendry et al. IL] based 
on coordinate transformation and Leonhardt [2[ based 
on conformal mapping method, invisibility cloak has at- 
tracted much attention and widely research. Notewor- 
thy, in 2003, Greenleaf et al. have suggested to de- 
sign anisotropic conductivities that cannot be detected 
by electrical impedance tomography through coordinate 
transformation of Poisson's equation which is similar to 
Pendry's method [3^,^. The effectiveness of transforma- 
tion based cloak with passive objects inside has been veri- 
fied through ray tracing approach [5, 6], full- wave simula- 
tions [7] , and analysis based on scattering models J8l-[ll| . 
On the other hand, the case that active sources are lo- 
cated inside the hidden area has also been investigated 
both in physics [12] and mathematical senses [13]. These 
investigations revealed the ideal invisibility cloak at par- 
ticular frequency can prevents the electromagnetic waves 
generated by the sources inside the cloak from propagat- 
ing out, meanwhile an extra surface electric and magnetic 
voltage are induced by the infinite polarization on the 
inner surface. As a result, it seems that the invisibility 
cloak with passive or active objects inside can not be de- 
tected by means of electromagnetic waves, while another 
method was suggested to detect the cloak by shooting a 
fast-moving charged particle through it [14]. 

In previous researches, the hidden passive objects and 
active sources are all electric neutral. As far as we know, 
the situation that the hidden objects are charged has 
never been discussed. However, the charges inside the 
cloak would be an important factor influencing the de- 
tection of the hidden objects. In this paper, we will only 
discuss the situation that the distribution of charges does 



* |zhangruoyang@gmail . com] 
"I" qzhaoyuping@bit.edu.cn] 



not change with time, so the shielding effect of the invisi- 
bility cloak reduces to an electrostatic problem. The elec- 
trostatic problem could be regarded as the limit of har- 
monically varying fields as cj ^ 0. However, some special 
property of harmonic fields could remain and the electric 
and magnetic fields would not decouple completely even 
at the limit cj ^ 0, such as \E / B\ = c does not varying 
with uj for a plane wave propagating in vacuum. So we 
treat the problem through solving the Poisson's equation 
strictly rather than taking the limit of zero frequency. 

In this paper, we fist construct the electrostatic cloak 
according to the invariance of Poisson's equation under 
coordinate transformation. Actually, the method to de- 
sign electrostatic cloak has been proposed by Greenleaf 
et al. [3, 4]. However, Greenleaf 's device is made of con- 
ductivity , while the cloak suggested here is composed 
of dielectric. Then, we calculate the fields in the whole 
space where the charges with arbitrary distribution are 
located inside a spherical cloak which are constructed by 
arbitrary radial transformation /(r). The result reveals 
the fields out of the cloak are only determined by the 
total charges inside the cloak and the shielding behavior 
of the cloak are very like the effect of charge confinement 
by a dual superconductor under the monopole-existing 
supposition, which model is often used in color confine- 
ment [15]. Moreover, for ideal cloaks, the result sug- 
gests an infinite electric energy when the total charges 
are nonzero. To overcome the embarrassment, we in- 
vestigate the unideal cloaks and find that they can also 
realize a good effect of shielding. 



II. ELECTROSTATIC CLOAK 

In vacuum, the Electrostatic fields obey Poisson's 
equation [sl, H- 
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where i/j is the electrostatic potential, 7*-^ is the con- 
travariant component of the spatial metric under an ar- 
bitrary curvilinear coordinate system (marked as S sys- 
tem) of the 3-D vacuous virtual space (VS), and 7 is the 
determinant of jij. In Eq. ([T]), y^7*-^ can be alterna- 
tively interpreted as the permittivity e'^^ of a real ma- 
terial expressed in another coordinate system of a real 
physical space (PS). Meanwhile, the physical quantities 
in VS also can be expressed in the same system as in 
PS, which is marked as S' system. Concerning an ar- 
bitrary radial transformation between S' and S system: 
r' = f{r)^ 0' = 6>, (j)' = (/), where S' is a spherical coor- 
dinate system, we obtain the relative permittivity of the 
spherical electrostatic cloak through the same process as 
in optical cloaking: 
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where we use e^^^^ to denote the uncoordinate compo- 
nents of the permittivity in orthonormal bases. 

We can see that the expression of permittivity of an 
electrostatic cloak is identical with it of an optical cloak 
[11]. So the electrostatic cloak can be regarded as an 
optical cloak at zero frequency. We will verify that the 
spherical electrostatic cloak expressed in Eq. (J2]) with "in- 
visibility conditions" /(a) = 0, f{h) = b can not be de- 
tected in electrostatic sense, where a, b are the inner and 
outer radiuses of the cloak respectively. 

Consider an electrostatic cloak wrapping a homoge- 
nous medium ei is located in vacuum endowed with an 
uniform electric field. A spherical coordinate system is 
constructed whose origin is at the center of the cloak and 
z axis is along the direction of the field. So the original 
field is i/jQ = —Eor cos 6> + C, where £^0 is the intensity of 
the uniform field. For charge-free anisotropic media, the 
potential satisfies V • ( £^ • V?A) = 0. Substituting Eq. (|2|) 



into the formula, we have 
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which is just the Poisson's equation in S' system of VS. 
Through separating variables, ip can be represented by 
series 
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where ij)^^^ ^ '^^^ ^/^^^^ represent the refiective field outside 
the cloak, the field in the cloak shell and the field in- 
side the internal hidden area respectively, and Pn{cosO) 
denotes the n-order Legendre function. 

According to the boundary conditions, which is the 
continuity of tp and the normal component of D across 
the inner and outer surface of the cloak, we obtain 
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and other coefficients are all zero, where 
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When /(a) = 0, f{b) = b, they reduce to a = —1, /? = 0. 
Therefore, ^°"* = 0, ^'"' = C, and the field in the cloak 
layer is 
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The result indicates the perfect cloak and its hidden 
medium can not be detected in that there is no reffec- 
tive field. Furthermore, the hidden area can not perceive 
the fields outside the cloak either. 

In terms of Eq. ([7j), we get the electric intensity and 



the electric displacement vector in the cloak layer: 
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The electric field intensity can be regarded as the tangen- 
tial vector of the electric lines of force, that is df/d\ = 
E'^^ where A is the parameter. Thus, we can derive the 
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FIG. 1. Electrostatic cloak in uniform electric field. The dis- 
tribution of potential are illustrated in y = plane. Red lines 
denote equipotential surfaces, (a) green lines denote electric 
lines of force, and (b) blue lines denote lines of electric dis- 
placement. The inner and outer radiuses of cloak are a = Im 
and 6 = 2m respectively. 



analytical expression of the electric lines of force 
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which employ the same expression of wave-normal rays 
in optical cloak exposed in a plane wave. Similarly, the 
lines of electric displacement are expressed by 



f{r) sinO = const, 



(10) 



It is just the expression of straight lines in S' system of 
VS, and it corresponds to the light-ray in optical cloak. 
The distribution of the fields are shown in Fig. [TJ where 
the constant of the potential are selected to be (7 = 0. 



III. ELECTROSTATIC SHIELDING 

Now we discuss the problem of electrostatic shielding of 
the cloak with an arbitrary distribution of charges inside. 
The simplest situation is that a point charge lies on the 

I 



z axis with a displacement ao(< a) from the origin. The 
potential in each region holds the same form as in Eqs. @ 
except that the total potential inside the hidden area 
should add the potential produced by the point charge 
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where q = q/{4:7rei) and q is the quantity of the point 
charge. The zero-point of potential is selected at infin- 
ity. When r > ao, the potential of point charge can be 
expanded as 
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Substitution the potentials into the boundary conditions 
at r = a and b yields 
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By solving the equations, all coefficients are determined 
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For the ideal situation /(a) = 0, /(6) = 6, the results are 
simplified to a„ = and /3o = £i/so, j3n = {n j^ 0). 
Thus the potentials reduce to 
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The field outside tiie cloak is precisely equal to the field 
produced by a point charge q at the origin in vacuum. 
While the potential in the whole of the hidden area tends 
to infinite since the term eiq/{sof{a)) exists as f{a) -^ 0. 
We can further solve the Green's function to this set 
of boundary conditions in light of the above solutions. 
The Green's function G(r, r ') is regarded as the potential 
produced by an unit point charge located at (r^, 0' ^ (j)'). 
Thus G(r, f') can be directly transformed from Eqs. (p!6|) 
through the substitution 
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For an arbitrary distribution p{f) of charges inside the 
cloak, the potential in the whole space is written as 



iPir) = J G{r,f')pif') 



r ') dV\ 
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where the domain of integration is the whole hidden area. 
The potential outside the hidden area still takes the ex- 
pression shown in Eq. (|16ap . which means the fields out- 
side is not affected by the distribution of charges inside 
the cloak but is only determined by the total charges 
q. Taking account of Gauss's theorem, it is actually not 
marvelous that the electrostatic cloak can not screen the 
outside space from the electric field of charged bodies in- 
side, as a result the charged bodies can be detected out- 
side the cloak. Nevertheless, apart from the total charges, 
no more information about the charge distribution can be 
detected, which is the least information allowed to gain 
under the restriction of Gauss's theorem. When the to- 
tal charges tend to zero, no field can go out of the cloak 
which is in complete agreement with the result for an 
active source inside the cloak as its radiation frequency 
goes to zero [12]. 

The behavior of electrostatic shielding of the cloak is 
identical with a spherical conducting shell in the screened 
area no matter where the sources are located, inside 
or outside the cloak. However, the electric field in the 
source-existing area is different between the two systems. 
It is interesting to generalize the two systems as two dif- 
ferent limitation of an unified system, a simple system 
with a dielectric sphere with permittivity ei placed in 
another medium 62- For simplicity, we also consider that 
a point charge q is located in the sphere with a distance ao 
from the center. The electric potential inside the sphere 
can be written as [iq 



V^" 



Y^r^ — 2aor cos 6 -\- a^ 

En + 1 a, 

n=0 



q- 



■£2 



£1 



(19) 






^r"P„(cos0). 



When £2 ^ CO, the result is the same as the case of 
conducting shell. On the other hand, in the limit £2 ^ 
0, the result is identical with the solution for the cloak 




FIG. 2. Electric field of two point charges with opposite quan- 
tity located at a distance a/2 from origin on positive and neg- 
ative semiaxis of z (a) in free space, (b) inside an ideal cloak, 
and (c) inside a spherical conducting shell. 



shown in Eq. (|16bp except the different constant term 
which in Eq. (|16b[) is q[si/{eof{a)) — 1/a] and which in 
the limit of Eq. (fT9|) is the n = term qeo/{£2Ci)^ while 
they both tend to infinity indeed. Actually, for spherical 
interface, only the radial component e(^rr) has been used 
in the continuity of the normal component of D. Thus 
Eq. (p!9|) with 82 = and the ideal cloak, whose radial 
component e(^rr) tends to zero as r ^ a, give the same 
result of field inside the hidden area. 

For the system of ideal cloak, the electric field in the 
hidden area is finite, despite the potential in the whole 
hidden area is infinite. While the field in the cloak shell 
tends to infinity as r ^ a. In addition, because the 
material of cloak is linear, the electrostatic energy of the 
system can be calculated by jl7| 



W = - I V^^"*(r Op(r dV\ 



(20) 



According to the infinite constant inside the cloak, the 
electrostatic energy of the system tends to infinity for a 
nonzero total charges even if the term of self-energy is not 
included. The result reveals it would cost an infinite work 
to put charged bodies into the cloak, which is impossible 
in practice. However, if the total charges are zero, the 
infinite constant of potential tends to zero. The trouble 
of infinite work would disappear either. 

Figure [2] shows an example of two point charges with 
opposite quantity located on positive and negative semi- 
axis of z with a same distance a/2 from the origin re- 
spectively. When the two point charges are inside the 
ideal cloak, all electric lines of force do not intersect with 
the inner surface of the cloak, while the equipotential 



surface are orthogonal to the inner surface (Fig. [2Jd). In 
contrast, as in the spherical conducting shell, all electric 
lines of force are orthogonal to the inner surface which 
itself becomes an equipotential surface (Fig. [2t). 



IV. SURFACE VOLTAGE AND EQUIVALENT 
SURFACE MAGNETIC CURRENT 



When checking the solution for an ideal cloak in 
Eqs. (p!6]) . one can find that the potential is discontin- 
uous across the inner interface of the cloak. For an 
arbitrary distribution of charges, the potential differ- 
ence across the interface r = a can be calculated by 
All; = j AG{ar,r')p{r)' dV, where AG(af,r') is the 
difference of Green's function across the interface: 
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This discontinuity of potential comes from the terms 
B^^fir)-^""^^^ in V^^^. For the ideal case /(a) -^ 0, the 
coefficients B^ -^ 0, so the terms of /(r)~*^"^+^) should be 
vanish in the cloak layer. While a meticulous calculation 
reveals that the limit B^f{a)~^'^~^^'^ is towards to a finite 
quantity which acts as a surface voltage to balance the 
potential difference Atp across the interface. The discon- 
tinuity also appears on the tangent component of electric 
fields. As we have pointed out in Fig.[2l the electric fields 
are tangent to the interface r = a on the inner side, yet 
are zero on the cloak side. Similarly, in the case of a 
cylindrical cloak with a transverse-electric (TE) incident 
wave, the tangent component of electric field is also dis- 
continuous at the inner interface of the cloak [10| . How- 
ever, the discontinuities are caused by different reasons 
in the two cases. In cylindrical cloak, the tangent com- 
ponent of magnetic field B(^^^ ^ oo at the inner surface 

r = a of cloak, and the integral ^ J^_q Bi^^j^^dr is equal 
to a finite value which acts as the surface magnetic dis- 
placement current to balance the difference of E^^^^ in the 

Maxwell's Eq. § E • dl = —- ^f B • ds for an infinitesimal 
contour across the interface |1QJ. While, for the ideal case 
of the shielding phenomenon, the normal component of 
electric field £^(^) becomes a delta function compressed 
on the interface and contributes to the surface voltage 

In the sense oiF dielectric materials, the surface voltage 
is caused by the infinite electric polarization of the mate- 
rial with e^^rr) = on the interface [12]. However the sur- 
face voltage and the shielding effect of charges can also be 
caused by the the surface magnetic current in the dual 
superconductor. As we all known, the Meissner effect 
says that a superconductor expels magnetic fields from 
its interior and therefore the equivalent permeability of 
superconductor is equal to zero. Similarly, in dual su- 
perconductor model, magnetic monopoles are supposed 
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FIG. 3. Distributions of (a) surface voltage A-^ and (b) sur- 
face magnetic current a (only component exists) on the 
inner surface of the cloak for the charge distribution shown in 
Fig. 2 with normalized q. 



to exist and the roles of electric and magnetic fields are 
interchanged. Moreover, the dual Meissner effect tries to 
expel electric fields out of the dual superconductors and 
gives a zero permittivity equivalently [15]. In this case, 
the discontinuity of the tangent electric fields across the 
surface of the dual superconductor comes from the sur- 
face magnetic current a, which satisfies the boundary 
condition 



E'""^ X Cr 



(22) 



As a result, the effect of electrostatic shielding of the 
ideal cloak is very like the charge confinement by the 
dual superconductors, so the inner surface of the cloak 
can be also interpreted as a dual superconductor layer. 
Although the medium of cloak is anisotropic, which is dif- 
ferent from the simple dual superconductor, we still can 
construct the ideal electrostatic cloak with dual super- 
conductive components as same as the way to construct 




FIG. 4. Electric field of two point charges with opposite quan- 
tity located at a distance a/2 from origin on z axis and x axis 
respectively (a) in free space, (b) inside an ideal cloak; and 
the distributions of (c) surface voltage Aip and (d) surface 
magnetic current a on the inner surface for the charge distri- 
bution with normalized q 




z[m] 



z[m] 



FIG. 5. Electric field of (a) a point charge located on pos- 
itive semiaxis of z with a distance a/2 from origin, (b) two 
equal and opposite point charges located as same as in Fig. [2] 
for an unideal cloak constructed by the linear transformation 
function with S = 0.2. 



magnetostatic cloak using superconductors suggested by 
B. Wood and J. Pendry [18]. Figure 3 shows the curves of 
the surface voltage and the equivalent surface magnetic 
current varying with for the charge distribution shown 
in Fig. 2. Since the distribution is symmetric about z 
axis, the voltage is also symmetric, and the surface cur- 
rent is along (j) direction. Figure 4 presents another ex- 
ample of charge distribution (Fig. 4a, b). Fig. 4c shows 
the voltage distributing on the inner surface, and Fig. 4d 
shows the surface magnetic current on the inner surface. 



V. UNIDEAL CLOAK 

In previous section, we have mentioned that the elec- 
tric energy tends to infinite for the case of nonzero 
charges inside an ideal cloak, and it is impossible in 
practice. To ease this embarrassment, we would take 
into account the case of unideal cloak which is also ob- 
tained from the radial transformation /(r). We still let 
/(r) satisfy f{b) = 6, however its zero point is not at 
r = a but at a' = (1 — (5)a, when S ^ the cloak ap- 
proaches to an ideal one. In this case, the potential inside 
the cloak no longer has the infinite constant term, and 
is continuous across the inner interface. A further cal- 
culation reveals that the unideal cloaks still have good 
property of electrostatic shielding both when the total 
charges is zero and nonzero. Figure [5] shows the case 
of unideal cloak designed by the linear transformation 
/(r) = b{r - a')/{h - a') with 5 = 0.2. For the case of 
a point charge located inside the cloak but not at the 
center, the fields outside the cloak still retain a highly 
spherical symmetry just as the fields generated by the 
point charge located at the center of the cloak, as shown 
in Fig. [5^. For the case of two equal and opposite charges 
located inside the cloak, the fields are mainly concen- 
trated inside the cloak, while the fields leaking out is 
extremely few, as shown in Fig. [SJd. 

To measure the effect of shielding accurately, we con- 
sider the electric flux propagating out of the cloak. Con- 
cerning the system of two equal and opposite charges in- 



^1-0; 




FIG. 6. (a) Linear transformation /(r) and (b) the corre- 
sponding radial component of permittivity e(^rr) varying with 
r in the cloak region under different value of 6 from to 1. 



side the cloak as shown in Fig. [2l3, Gauss' theorem tells 
the total flux out of the cloak is zero. However, if we 
calculate the fiux out of two hemispherical surfaces of 
a sphere r = ro > b split hy z = plane respectively, 
the pair of flux must have an equal magnitude |^| but 
opposite sign because of the symmetry of charge distri- 
bution. Thus the absolute value of electric fiux out of 
each hemisphere would be a suitable quantity to mea- 
sure the effect of shielding for this system. We still con- 
sider the linear radial transformation /(r) with different 
S (Fig.[6ti) and corresponding permittivities (Fig.[6]3). In 
this situation, |^|/^o varying with 6 is shown in Fig. 7, 
where | ^ | is the absolute value of fiux flowing out of the 
hemisphere and ^o is the absolute value of fiux out of the 
hemisphere with no cloak existing. The slope of |^|/^o 
varying with S is equal to zero as (5 ^ 0. The result man- 
ifests that a nearly ideal cloak would still realize a good 
effect of shielding. In fact, the derivative of e^^rr) of ideal 
cloak with respect to r is equal to zero for arbitrary /(r), 
therefore a nearly ideal cloak with small value of 6 would 
always have a good behavior of shielding. The effect of 
shielding can be measured from another aspect of energy 
density ^E ■ D^ as shown in Fig. [51 We can see that the 
energy density outside a nearly ideal cloak tends to zero 
when 5 is very small and is much smaller than the case 
of no cloak existing. Along with the increase of 6, the en- 
ergy density increases nearly stable and goes to the limit 




FIG. 7. The absolute value of electric flux |^| divided by ^o 
flowing out of the hemispherical surface split by z = plane 
varying with S. 
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FIG. 8. Electric field energy density varying on z axis under 
different value of ^, where the red line denotes the case of no 
cloak existing (^ = 1). 



of no cloak existing when 5^1. To sum up, the slight 
breaking of the perfection would not change the property 
of electrostatic shielding for the cloak significantly. 



VI. CONCLUSION 

To summarize, we have verified that electrostatic 
spherical cloak takes the same form of permittivity as the 
invisibility cloak under nonzero frequency, and demon- 



strated its behavior of electrostatic shielding is identical 
with a spherical conducting shell in the screened region. 
If the electrostatic sources are outside the cloak, the field 
can not propagate into the screened region inside the 
cloak, on the other hand, if the charges with arbitrary dis- 
tribution are inside the cloak, the field outside the cloak 
is just as the field generated by a point charge located 
at the center of the cloak, so the only information which 
can be detected is the total charges inside the cloak. For 
ideal case, the potential across the inner interface of the 
cloak is not continuous in that £{^rr) — causes the infi- 
nite polarization on the inner interface. However, £ = 
can be also interpreted as the property of a dual super- 
conductor, and the behavior of shielding are also very 
like the property of charge confinement caused by dual 
Meissner effect. Another problem existing in ideal case is 
the infinite field energy when the total charges inside the 
cloak is not zero. Nevertheless, the problem no longer 
exists for unideal case, in addition the nearly ideal cloak 
also have a good effect of electrostatic shielding. 
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